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We analyze the stability of a gyroscopic oscillator interacting with a finite- and infinite- 
dimensional heat bath in both the classical and quantum cases. We consider a finite gyroscopic 
oscillator model of a particle in a magnetic field and examine the stability before and after 
coupling to a heat bath. It is shown that if the oscillator is gyroscopically stable, coupling to a 
sufficiently massive heat bath induces instability. The meaning of these ideas in the quantum 
context is discussed. The model extends the exact diagonalization analysis of an oscillator and 
field of Ford, Lewis, and O'Connell to the gyroscopic setting. 

1 Introduction 

In this paper we investigate the stability of a gyroscopically stabilized system interacting with a 
finite dimensional heat bath. More details will be given in a forthcoming paper [?]. 

A gyroscopically stabilized system is one that is unstable without gyroscopic forces but becomes 
stable with the addition of these forces. Infinitesimal dissipative perturbations are known to induce 
instability in Hamiltonian systems that are gyroscopically stabilized; see [?]. 

Since the origins of dissipation {e.g. friction, viscosity...) lie in the transfer of energy from 
one form (energy of one subsystem) to another form (that of a second subsystem) of a larger 
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Abstract 



conservative system, it is natural to expect the analogue of the above destabilization phenomenon 
to be present within the more fundamental context of conservative systems which exhibit internal 
energy transfer. In [?] and [?], we explore this in the context of a gyroscopic oscillating mechanical 
system coupled to an extended wave system (infinite string). Due to the coupling, motion within 
the mechanical system generates waves which can be carried off to infinity. Such radiation damping 
has been studied in models arising in the theory of quantum resonances, ionization type problems 
and nonlinear waves; for more detail see [?], [?], [?], [?], and references therein. 

As a model of a gyroscopically stabilized system, we consider a charged particle stabilized by 
a magnetic field. We show that even for a finite heat bath sufficient coupling strength induces 
instability. This result is quite striking in the sense that one does not expect a finite heat bath to 
mimic dissipation. A graphical criterion is given for determining the onset of instability. Further 
we generalize the analysis to the quantum setting. The infinite limit and oscillator susceptibility 
are analyzed in the forthcoming paper [?] mentioned above. 

In the quantum setting we show that while a stable oscillator has positive energy bound states, 
a gyroscopically stable oscillator has both positive and negative energy bound states. Coupling 
to the bath with sufficient coupling strength induces unbound states in the gyroscopically stable 
system. Applications of such quantum systems, such as Penning traps, are used to trap small 
molecules and obtain extremely precise measurements of atomic quantities [?]. 

Our analysis extends the heat bath analysis of [?] to the gyroscopic setting. We show here that 
there is a beautiful extension of their graphical (intersection-theoretic) criterion for stability to a 
more complex class of curves. 

2 Example of a Chetaev system 

In this section, we discuss a physical realization of the gyroscopic oscillators discussed in this paper. 
Consider the motion of a charged spherical pendulum in a magnetic field whose linearization is that 
of a charged planar oscillator in a magnetic field (for details of the full nonlinear system see e.g. 

[?]•) 

Let B be a divergence-free vector field. Let A be the vector potential, B = V x A. Note if 
we choose B to be the constant magnetic field in the direction normal to the plane of oscillation, 
the vector potential can be chosen as , A = ^B x q, where q = (x,y,0) T is the position of the 
oscillator. 

Assume the oscillator has unit mass and unit charge and that the speed of light is unity. The 
Lagrangian, L : TM 2 — > R, is defined by: 

£(q,q) = J||q|| 2 + A- q -[/(q) 

\ (2-1) 
U(c l ) = -(ax 2 +(3y 2 ). 

Choosing B to be of constant strength B, normal to the plane of oscillation, we obtain the Euler- 
Lagrange equations of motion: 



x — By + ax = 
y + Bx + f3y = 0. 



(2.2) 



Remark 1: If a and (5 are both negative, the oscillator in a field system is unstable for small B. 
However, if B 2 + a + (3 > 2y/a[3, the oscillator stabilizes, i.e. the eigenvalues are on the imaginary 
axis - this is what is referred to as gyroscopic stabilization. For further details see below and [?]. 

Remark 2: If the system is gyroscopically stable it can be shown that adding a small amount 
of dissipation to the system renders it unstable (i.e. there are unstable eigenvalues). For a more 
precise statement and generalizations see section ||. 

3 Stability of gyroscopic systems 

We recall here some general properties of linear systems with gyroscopic forces. The systems above 
are examples of such systems. 

The general form of a gyroscopic system is 



where q G M n , M is a positive-definite symmetric nx n matrix, S is skew, and V is symmetric. As 
in [?] we shall call this the Chetaev system (see [?]). 

We say the system is gyroscopically stable if for S = the origin is an unstable equilibrium, 
but for 5/0, the origin is a spectrally stable equilibrium (i.e. the eigenvalues of the linearized 
system have non-positive real part). The matrix S is sometimes referred to as a magnetic term 
which arises from charged oscillators in a magnetic field. 

An important property of this system is that it is normal form for a simple mechanical system 
about a relative equilibrium which is given modulo an abelian group. That is, it is the normal form 
of a system defined on the cotangent bundle TQ of a configuration space Q where the Lagrangian is 
given by kinetic minus potential energy. One can obtain a similar normal form in the case of a non- 
abelian group. (See [?] and [?].) The magnetic term naturally arises in the symplectic form when 
investigating the quotient space. However, we can obtain the same dynamics from a canonical 
symplectic form and an augmented Hamiltonian. This can always be done by the momentum 
shifting lemma (see [?]). Classically, the two representations of the Chetaev systems are equally 
useful. However, the canonical bracket is preferable when quantizing the mechanical system. (See 
[?], [?], [?] for further physical discussions.) 

As for gyroscopic stability, the number of negative eigenvalues of the quadratic form plays a 
crucial part as [?] discusses and which we summarize in the following proposition: 

Proposition 3.1 Consider the canonical gyroscopic system Mq + Sq + Vq = 0, where M is a 
symmetric positive definite matrix, S is a skew-symmetric matrix, and V is a symmetric matrix: 

• If V has a odd number of negative eigenvalues (counting multiplicity) then the origin is an 
unstable equilibrium. 

• IfV has an even number of negative eigenvalues (counting multiplicity), we can choose S so 
that the origin is a spectrally stable equilibrium. 

We omit the proof here. 

Gyroscopically stable systems exhibit interesting instability when perturbed by dissipative 
forces. Suppose now that V has at least one negative eigenvalue. A key result of [?] is that 
adding small dissipation always yields instability. More precisely it is shown that 



Mq + Sq + Vq = 




Theorem 3.2 Under the above conditions, if we modify the general Chetaev system by adding a 
small Rayleigh dissipation term, 

Mq + (S + eR)q + Vq = (3.2) 
for small e > 0, where R is symmetric and positive definite, then the perturbed linearized equations 

z = L e z, 

where z = (q,p) are spectrally unstable, i.e., at least one pair of eigenvalues of L e is in the right 
half plane. 

This result builds on basic work of [?], [?], and [?]. We refer to this as dissipation induced 
instability. We see that the Hamiltonian of a gyroscopically stabilized Chetaev system is indefinite. 
In this case, Rayleigh dissipation decreases the value of the Hamiltonian, but this does not bound 
the motion of the Chetaev system. In particular, the Hamiltonian may decrease to zero, while the 
displacement and velocities grow exponentially. 

[?] also prove a similar stability result for the non-abelian case, but the abelian result is sufficient 
for our purposes. 

4 Oscillator coupled to a bath 

We now consider the gyroscopic oscillator coupled to a bath of oscillators via an augmented La- 
grangian 

Lcoupled = L + Lbath (4-1) 

where L is defined in equation ( |2.1] ) and 

N 1 

Lbath = ^2 m i( i: ' 2 j +V 2 j) 

3=1 

--mjJI{{xj - xf + (yj - y) 2 ), 

with < L0\ < ■ ■ ■ < lon being the characteristic frequencies of the bath. Similarly, the Hamiltonian 
can be augmented in corresponding fashion. 

This model extends the model described by [?] to the gyroscopic setting. As discussed in that 
paper this model provides a good physical realization of oscillator-bath coupling. For simplicity of 
the computations here we restrict ourselves to the case a = (3. Similar results hold in the general 
case. We also restrict ourselves to the generic situation though again similar results hold in the 
non-generic setting. 

The theorem below in fact generalizes the results of [?] in a rather beautiful fashion, leading to 
the study of a more complex intersection problem. 
We can show: 



Theorem 4.1 Consider the case where a = f3 and the equation 

In 1 o w 2 w 2 

( W ± -5) 2 - (a + = £ m,^^ (4.2) 

3=1 J 

Generically we have the following: 

The oscillator is stable if there exists AN + 4 rea/ frequencies u = Qi ,i = ±1, • • • db (2iV + 2) 



which are solutions to the equation (J^.i). Let VL\ be the smallest positive solution and let be 
increasing with respect to the index i. 
More precisely we have 

(i) In the case of instability of the oscillator, (a + \B 2 < 0), there are only AN real solutions. 
Additionally, there are 2 pairs of imaginary solutions corresponding to instability. 

(ii) In the case of strong stability (a > 0), there are AN + A real, normal modes maintaining 
stability via coupling. 

(in) In the case of gyroscopic stability (a + \B 2 > 0), there are two possibilities. The values of 
rrij and ujj determine which occurs: 

(a) If Oi < u\, then there are AN + 4 real, normal modes (counting multiplicities) corre- 
sponding to stability. 

(b) If ill > ^l) then there are AN real, normal modes and 2 pairs of imaginary mode 
corresponding instability in the coupling. This case occurs for large rrij . 

Idea of Proof. The equations of motion are 

TV 

x — By + ax = mjLo'j(xj — x) 

3=1 
N 

y + Bx + (iy = ^ m j u;'j(y j - y) 

3=1 

Xj + UJjXj = Uj'jX 

Vj+UjVj = u)y. 
We seek normal mode solutions of the form 

x{t) = X Q {u>y ut ; y(t) = Y Q {u>)e iut 

x j (t)=X j (u)e i » t ; y j (t)=Y j (u)e i " t . 

In the isotropic case (a = f3), standard normal mode analysis yields nontrivial solutions for the 
characteristic equation ( |4.2j ). Stability can be insured if all normal modes are real. The precise 
results follow simply from examining the number of real intersections of the curves on the left and 



right sides of equation (f4.2j). For all cases it is helpful to refer to figures 4.1, 4.2 for the case N = 1 
plotting left hand side and right hand side of equation (|47 

Details of the proof are given in [?] . 
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Figure 4.1: Stable, a 

5 Gyroscopic quantum oscillators 



= —1, B = 3, mi = j, and uj\ = 1. 



The instability results for the classical gyroscopic oscillator coupled to a heat bath can be extended 
directly to the quantum setting. This follows from the general quantization procedure for a classical 
system with quadratic Hamiltonian and gyroscopic terms. We give firstly the general analysis here 
as this is useful for understanding how the diagonalization procedure should be applied in the 
quantum case. 

Consider a general classical Hamiltonian of the form 



H(q, p) = \p T M- X p + q T Sp + -q T Kq 



(5.1) 



where q,p are the n-dimensional position and momentum vectors respectively, and M,S,K are 
constant n x n matrices with M, K symmetric and S skew-symmetric. For 



J 



/ 
-/ 



where / is the n x n identity matrix, the Hamiltonian equations of motion are given by 

= JVH = 



dt 



-Sq + M - V 
— Sp — Kq 



(5.2) 



(5.3) 



Now consider the quantization of the above system. The Heisenberg equations of motion are given 
by 



ihq = [q, H] 
ihp = [p, H] 

where qi,Pi obey the standard commutation relations. This gives 

-Sq + M~ x p 
—Sp — Kq 



in— 
dt 



q 




P. 


= i 



(5.4) 



(5.5) 



Figure 4.2: Unstable, a = — 1, B = 3, mi = 2, and wi = 1. 



just as in the classical setting as an easy computation shows. 

The point is that the eigenvalue computation for the equations of motion in the classical setting 
is precisely equivalent to that for the Heisenberg equations of motion. Thus for the oscillator 
coupled to the heat bath the numerical computation in the classical case gives us the correct 
eigenvalue information for the quantum case also and we can deduce qualitative behavior as for the 
single gyroscopic oscillator in the previous section. Thus, for sufficiently large coupling, we observe 
unbound states even when the uncoupled gyroscopic oscillator exhibits bound states. 

One also observes that diagonalizing the Heisenberg equations of motion is in fact equivalent 
to diagonalizing the Hamiltonian. 

2-D gyroscopic quantum oscillator 

The quantization of the uncoupled gyroscopic system illustrates that the stability analysis of the 
quantum system is analogous to the classical normal mode calculation. However there are some 
subtleties in the interpretation of the stability analysis and for this reason we consider in detail 
the spectral analysis of an uncoupled quantum gyroscopic systems. By the analysis of the previous 
section similar considerations apply to the full system of the oscillator coupled to the heat bath. 
We follow the Dirac formalism as described in [?] for example. 

Analogously to the classical coupled oscillators, we define the quantum Hamiltonian H for an 
oscillator in a magnetic field of strength B and isotropic potential, 



H 




+ ^(PxV ~ PyX). 



For a quantized gyroscopically stable oscillator, we use a standard change of coordinates to 



write the Hamiltonian in terms of normal modes: 




where ujq = \J\^j- + ot\ and properties that [j4 r ,A s ] = [ylJ,A|] = and [A r ,Aj] = 6 rs for r,s £ 
{+, — }. In the new coordinates, we have 

H = % -f)(4^+ + ^) + 
^( Wo + |)(^lA_ + i). 

We verify that the new frequencies ujq ± ^ are exactly the classical frequencies from the normal 
mode calculation. In the Heisenberg picture, we have 

ihAl = [Al,H] = -h(uj T^)Ai 
A\_ = e i{u > oT f )t Al(t = 0). 
The eigenvalues of the Hamiltonian, E ni>rl2 are 

E num = hwQ{m + n 2 + 1) + h—(m - n 2 ), 

where ni, n 2 € {0, 1, 2, • • • }. In the case of gyroscopic stability, we have > ^0; and the energy 
spectrum is indefinite and unbounded above and below. In the case of pure stability (a > 0), we 
have -j < ojq, and the energy spectrum in positive definite. 

In the case of a quantized unstable gyroscopic oscillator, (y + a < 0), it can be shown there 
are no bound eigenstates. 

Similarly in the case of a quantum oscillator coupled to a finite heat bath we can use the change 
of variables in equation ( p^q) to compute the characteristic frequencies, thus reducing the coupled 
problem to a system of independent quantum oscillators. Thus the spectrum of the Hamiltonian 
of the coupled quantum system is positive if and only if the classical system is strongly stable; the 
spectrum is discrete with both positive and negative eigenvalues if and only if the classical system 
is gyroscopically stable; and finally, we have unbounded eigenfunctions if and only if the classical 
system in unstable. 



6 Conclusions 

Often one models dissipation or friction by coupling a mechanical system to a thermal reservoir. We 
have modeled the thermal reservoir as a collection of oscillators coupled to the mechanical system. 
If the reservoir is sufficiently massive, then the coupling may induce instability in the mechanical 
system. Equivalently, if the reservoir contains sufficiently low frequencies, then similar instabilities 
may arise. 



Unlike other models of dissipation (i.e. Rayleigh), the conservative nature of the thermal reser- 
voir allows us to naturally extend our results to the quantum setting. We interpret stability in 
the quantum setting as corresponding to bound eigenstates. A gyroscopically, stabilized, quantum 
Chetaev system exhibits this stability while having an unbounded spectrum. 



